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We investigate dynamic hysteresis in ferromagnetic thin films with zigzag domain walls. We 
introduce a discrete model describing the motion of a wall in a disordered ferromagnet with in-plane 
magnetization, driven by an external magnetic field, considering the effects of dipolar interactions 
and anisotropy. We analyze the effects of external field frequency and temperature on the coercive 
field by Monte Carlo simulations, and find a good agreement with the experimental data reported 
in literature for Fe/GaAs films. This implies that dynamic hysteresis in this case can be explained 
by a single propagating domain wall model without invoking domain nucleation. 
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I. INTRODUCTION 

Ferromagnetic materials are concrete examples of co- 
operatively interacting many-body systems. When a 
magnet is driven by a varying external magnetic field, the 
system cannot reach equilibrium instantaneously, due to 
the internal relaxation delay. If the external field oscil- 
lates, the magnetization will do so as well, lagging behind 
the field. This effect gives rise to a non-vanishing area 
of magnetization-field loop whose form will depend on 
the applied field frequency: a phenomenon known as dy- 
namic hysteresis^'^. Since the loop area represents the 
amount of externally supplied energy that is irreversibly 
transformed into heat during one magnetization cycle, 
dynamic hysteresis has important technological implica- 
tions, e.g. for high frequency devices applications. Fur- 
thermore, from a purely theoretical point of view, the 
dynamics of disordered magnetic systems represents a 
central problem in non-equilibrium statistical mechanics. 
While in three dimensional systems dynamic hysteresis 
is well understood in terms of eddy currents dissipation, 
this effect is expected to become negligible by reducing 
the sample thickness'^. Thus recently there has been a 
renewed interest in understanding two dimensional sys- 
tems, both experimentall3r2i4i5,6,7,8,9,io,ii,i2,i3,i4,i5,i6,i7 

and theoreticalljii*i2ii2i2£, motivated by the applications 
of thin ferromagnetic films in magnetic recording tech- 
nology and spintronic devices. 

Two classes of models are mostly used to investigate 
the magnetization reversal properties on a microscopic 
scale, spin models of the Ising typ9i*22iSi24, or extended 
domain wall modelsSSiSiiS^. The theoretical tools used to 
interpret experimental data on dynamic hysteresis are of- 
ten grown out of the first class of models which suggest a 
universal scaling law for the dependence of the hysteresis 
loop area A on the external parameters, i.e. the tem- 
perature of the system T and the amplitude and fre- 
quency uj of the applied magnetic field. In particular, it 
is expected from the models that A oc HqUJ^T~'^ , where 
a, P and 7 are the scaling exponents^. The experimen- 
tal estimates of these exponents display, however, a huge 
variability^i'^i^i'^i^i^d°i"d^d^ and the validity of that uni- 



versal scaling law is still under debateiStS^. Since various 
phenomena may in principle contribute to the hysteretic 
behavior, like domain nucleation, domain wall propaga- 
tion or simply retardation of the magnetization due to 
fluctuations, it has been proposed that the lack of good 
scaling of the function A{u!) is due to a cross-over be- 
tween two distinct dynamical regimes, one dominated by 
domain wall propagation, and the other by nucleation of 
new domains^ 

The second class of models used to investigate ferro- 
magnetic systems considers the dynamics of individual 
domain walls as the relevant mechanism for hysteresis. 
In two dimensional systems, developing such a kind of 
model can be even more complicated than in the bulk 
three dimensional case, due the possibility for the mag- 
netization to lay in or out the film plane, and the huge va- 
riety of domains and domain walls topologies (for an ex- 
haustive overview of the existing configurations together 
with many experimental images, see ref. ^2^. Dynamic 
hysteresis due to the motion of 180° Bloch domain walls 
has been extensively investigatediSiSS, but less is know 
about charged walls. 

In this article we focus on two dimensional sys- 
tems with zigzag domain walls, arising from the com- 
petition between dipolar forces and magnetocrystalline 
anisotropy in thin films with head-on magnetization be- 
tween nearest-neighbor domainsS^. These walls have 
been originally observed in thin film magnetic record- 
ing media, where head-on domains are induced by means 
of the application of a recording head field, and have 
been then reported in several magnetic materials such 
as ironSS, Co2£, Gd-Co2&, epitaxial Fe films grown on 
GaAs(OOI)^, finemets and many others. In addition, 
zigzag walls have also been observed in ferroelectric mate- 
rials, such as Gd2(Mo04)3 crystals'^^. Most calculations 
reported in literature for zigzag walls focus on the deriva- 
tion of the equilibrium parameters (e.g. zigzag angle and 
amplitude or period)^^*^ and do not consider their dy- 
namics. 

Here we introduce a discrete model for the motion of 
a single zigzag wall in a disordered ferromagnetic two 
dimensional sample with in-plane uniaxial magnetiza- 
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tion, driven by an external (triangular) magnetic field. 
The model is based on the interplay between dipolar 
and anisotropy energy contributions, in the presence of 
structural disorder. Dynamic hysteresis is investigated 
by Montecarlo simulations analyzing the behavior of the 
coercive field He as a function of the external field fre- 
quency, temperature T, and sample thickness. We find 
good quantitative agreement with experimental data re- 
ported for Fe/GaAs thin films^. Our results indicate 
that the experiments can be interpreted by a domain 
wall propagation model, and thus ruling out explana- 
tions involving a cross-over with domain nuclcation or 
other processesSiS. 

The manuscript is organized as follows: in Sec. ITU we 
present an overview on the energetics of a zigzag domain 
wall, computing magnetostatic, anisotropy and disorder 
energies. In Sec. IIIII we estimate the mean zigzag half- 
period and the coercive field and compare our result with 
experimental observations, in order to test the reliability 
of our approximations. Next, in Sec. IIVI we present our 
model and the results obtained by Montecarlo simula- 
tions for the frequency, temperature and thickness de- 
pendence of the coercive field, and compare them with 
experiments. Our results are finally resumed in Sec. IVII 



II. ENERGETICS OF ZIGZAG DOMAIN WALLS 

In thin uniaxial ferromagnetic films, we can distin- 
guish between two main classes of domain walls: the first 
is represented by prevalently straight (magnetically un- 
charged) walls parallel to the easy axis, and the second by 
charged walls separating two domains with head-on mag- 
netization. Since there is a cost of magnetostatic energy 
associated with the magnetic charge that increases with 
the sample thickness, these walls are observed mostly in 
thin filmSi. On the other hand, a charged straight wall 
is unstable and becomes metastable by forming a zigzag 
pattern to minimize its energy. 




FIG. 1: The magnetic configuration in the Neel tail. 

A quite extensive derivation of the equilibrium zigzag 
parameters (amplitude, period and angle) has been re- 
ported in Ref. y^- The calculation is based on a Neel 
tail transition model^i which involves an in-plane mag- 
netization rotation over a transition region: the spins 
rotate following the wall shape. Thus the entire region 



surrounding the wall exhibits a nonuniform magnetiza- 
tion (see Fig^. This spreading occurs at the cost of 
increasing the anisotropy energy. The total energy may 
be expressed as a function of the zigzag parameters that 
are then obtained by minimization. From our point of 
view, it is important to stress that from magneto-optical 
images it can be inferred that the zigzag angle is constant 
across the wall and does not change during the motion. 
Thus, in the following, we will consider the angle as a 
fixed parameter in our model. As we mentioned in the 
introduction, the zigzag shape of the wall is due to the 
interplay between the magnetostatic and the anisotropy 
contributions to the total energjE-^^. The magnetostatic 
term opposes a straight wall, which would maximize the 
magnetic charge density, and favors a large zigzag ampli- 
tude, so that the magnetic charges at the wall (all of the 
same sign) are as separated as possible. The anisotropy 
term prevents the amplitude to increase freely, avoiding 
an excessive deviation of the magnetization from the easy 
axis associated with a spread out Neel tail. 

Our purpose is to study domain wall motion under an 
external magnetic field, by discrete model simulations. 
To this end, we calculate the total energy of an arbitrary 
zigzag wall configuration. As we are interested in the 
macroscopic response, we do not consider the details of 
the wall internal structure, and treat only the magneto- 
static, the anisotropy and the disorder contributions: 

E = Em + Ean + Edis- (1) 

In equation J^l, the magnetostatic term Em takes into 
account the interaction between magnetization and stray 
field, the anisotropy Ean is the energy cost of deviations 
from easy axis and Edis models structural disorder, im- 
purities, defects and so on. In the following subsections 
we will discuss these terms in more detail. 



A. Magnetostatic energy calculation 

We consider two generic segments of a zigzag wall of 
total length L. We label the segments as i = l,...,rt, 
where n = L/p and p is the half-period of the zigzag. We 
call h the zigzag amplitude and 6 the angle between the 
zigzag segment and the easy axis. The thickness of the 
film is e, and it coincides with the wall thickness as we 
will consider only rigid walls (see Fig. |21for a definition 
of the parameters). 

Since the magnetostatic self-energy only depends on 
the total magnetic charge, which is constant during the 
wall motion, we limit the calculation to the magnetostatic 
interaction energy (i.e. i ^ j). The contribution due to 
the i-th and the j-th segments may be written as the 
surface integral 

Eij ^ j a<j)^{x,y,z)dSj , 

where cr = M • n = 2MsSin0 is the (constant) surface 
magnetic charge density {Ms is the saturation magne- 
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FIG. 2: (Color online) Sketch of the parameters of the zigzag 
wall. The easy axis is along the y direction. 



tization), n is the unit vector normal to wall segment 
surface and Sj is the surface of the j-th segment. The 
scalar potential generated by the i-th segnrent, 4>i{x, y, z), 
is calculated by 



(j)i{x,y,z) = 




So, for e small with respect to the segments distance, we 
can write 

E,j = SM^eVo / dx / dx' 

J jp J ip 

1 

— x')^ + {rriiX + qi — rrijx' — qj)^ ' 

(2) 

where mi and nij are the slopes (which values should be 
ih/p = ±tan(7r/2 — 9) ) and qi and qj the j/-intercepts 
of the i-th and j-th segments. The direct solution of 
equation ij^Jl is very involved, and we report it in the 
appendix. 



B. Anisotropy energy calculation 

The anisotropy energy term Ean describes the energy 
cost of the deviation of magnetic moments from the easy 
axis of the material, which in the simple case of an uni- 
axial crystal can be written as 

Ean = Jd^rKu sin^ (j), (3) 

where is the in-plane uniaxial anisotropy constant 
and (f> is the angle between the easy axis and the mag- 
netization vector. The rotation of the magnetization is 
associated with the Neel tail (FigQ . We assume^^ that 
the charge is uniformly distributed over the entire band 
containing the wall. Although this approximation exag- 
gerates the diffuseness of the charge, it has been used 
to calculate vertex angles and zigzag amplitudes that re- 
sulted to be in reasonable agreement with experimental 
observations^. We can estimate Eq. O (with the nota- 



tion sketched in Fig. [^J and obtain: 

Ean = , . / dy I dxsiv? (j){x,y). 

fttant/J_,j/2 J -(h/2-y)ta.n9 

(4) 

Assuming that (l){x^ y) describes a linear in-plane 
rotation of the magnetization vector, (/)(a:, y) ~ 
9x/ [{h/2 — y)ta,'a.S\^ developing sit? (j) in power series 
and then integrating term by term, we obtain the 
anisotropy energy for unit length as 

Ean = eKuhc{B), (5) 

where c{9) is a constant function of the zigzag angle 9: 

(^_iynj Q2(m+l + l) 

- (2m-M)!(2; + 1)! 2im + l + l) + l' 

ni,l—0 ^ 

which could be evaluated numerically. 



C. Disorder 

Different sources of inhomogeneities are found in vir- 
tually all ferromagnetic materials, and the presence of 
structural disorder is essential to understand the hys- 
teretic behavior, especially to account for the residual 
coercive field at w — > 0. Disorder is provided by va- 
cancies, non-magnetic impurities, dislocations or grain 
boundaries in crystalline systems, variations of the easy 
axis between different grains for polycrystalline materi- 
als and internal stresses for amorphous alloys. We will 
consider only quenched (frozen) disorder, that does not 
evolve on the timescale of the magnetization reversal. For 
simplicity we model disorder by an energy contribution 
associated to each site which may be occupied by a seg- 
ment (our discrete unit length) of the zigzag wall. This 
term is extracted from an uncorrelated random Gaussian 
distribution with zero mean. 

The dependence of the disorder energy on the film 
thickness can be obtained by simple considerations. 
Given that the minimal segment has length a, we dis- 
cretize the sample thickness in e/a elements and associate 
a Gaussian random variable to each of them. The mean 
square value of the disorder energy per unit domain wall 
length is then proportional to the sum of contribution 
over the whole thickness, and thus to e. As a consequence 
of this, the disorder contributes to the total energy (eq. 

by a term proportional to the square root of the film 
thickness 

Ed,s = V^A, (6) 
where A is the variance of the random variables. 
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III. THEORETICAL CONSIDERATIONS AND 
COMPARISON WITH MATERIALS 

A. Zigzag parameters 

A way to test the reliability of our model is to compare 
the results of our model for some relevant parameters 
with the measured experimental values. As an exam- 
ple, we can estimate the typical zigzag half-period Peq for 
an equilibrium configuration in absence of external field. 
An approximation of the magnetostatic energy can be 
obtained in closed form by developing Eq. [21 for p « L 
and is given, for unit length, by 



En 



In the absence of disorder and for Hext = 0, the total 
energy (Eq. ^ can be written as 



E = E„, + Ean ^ SM^e^fio In (L/p) + eK^ 



P 



tan 6 



where we have imposed h = p/ taxiO. If the configuration 
is stable, considering the T ~ case, we can impose 
dE/dp = and obtain 



tan ( 



= 0, 



from which follows 



SM^efio tan 



(7) 



We can estimate the numerical value of p^g by using 
the parameters reported in literature for typical ferro- 
magnetic thin films. For example, for Fe/GaAs{00l) 
analyzed in Ref. we can set /igMg — 2T and — 
0.5 X 10^ J/m'^, so that for a thickness e — 25nm and an 
angle 9 = 20°, we obtain 



Pe 



100 /xm, 



which is in good agreement with the typical lengthscale 
inferred from magneto-optical investigations^. 



B. Coercive field and thiclcness dependence 

Another quantity that is interesting to obtain quanti- 
tatively is the value of the (zero frequency) coercive field. 
A very rough order of magnitude estimate could be ob- 
tained supposing that the disorder is small and that the 
anisotropy and magnetostatic terms are of the same or- 
der of magnitude. Close to the coercive field when the 
energy variation is stationary we can set 



HoMsephHc = ephKuc{0), 



which implies 



(8) 



Employing the parameter values reported above, we ob- 
tain He ~ 15 Oe, which is in reasonable agreement with 
the range of values found for Fe/GaAs thin filma^. 

A more refined quantitative calculation of He is diffi- 
cult because we can not easily quantify the value of the 
disorder term. Nevertheless, we can obtain from simple 
considerations the thickness dependence of coercive field. 
The energy difference between two zigzag wall configura- 
tions in the presence of external magnetic field H^xt is 
given by (see eq. 

AE ^ AE,n + AEan+ AEdis+AEext, (9) 

where the interaction AE^xt with H^xt is given by 



AE, 



-2fioM,eHextAA, 



(10) 



and is the area interested by the magnetization re- 
versal. As it could be seen from equations 130 and El eq. 
could be rewritten as 



AE- 



e^AEL 



eAE^,, + VeAE',,, 



(11) 



where AE',^ = AE„,/e'', AS^„ = AE^n/e, 
AEdisI Ve and AE'^^^ — AE^xt/^ do not dependent on e. 
As AE represents the energy barrier that the zigzag wall 
has to overcome in order to move in the direction of the 
external magnetic field, we expect that Eq. encodes 
crucial informations on the dynamics. 

The coercive field dependence on the thickness, can 
be obtained comparing AEm with the most relevant of 
the various terms contributing to AE. For e 0, 
AE is dominated by the disorder contribution. Com- 
paring AE„i (X with AE^is oc ^/e, we easily obtain 
He oc 1/ \/e. If the wall is strongly pinned by disorder or 
anisotropy, starting from a saturated configuration, small 
external field changes are not able to trigger avalanches, 
resulting in square-shaped hysteresis loops with high co- 
ercitivity. Otherwise, if the pinning is weak, avalanches 
are induced even by relatively small field changes, so that 
the loops will be tighter, with small coercitivity. We can 
expect that the 1 / -^/e dependence will disappear in the 
limit of weak disorder regime. For intermediate e values, 
the leading terms in eq. mwill be the anisotropy energy 
AEan OC e and the interaction with the external magnetic 
field AEexti leading to a thickness independent coercive 
field given by Eq. [HI Finally, for larger values of e, the 
dipolar energy will lead the energy barrier balance, in- 
ducing a linear dependence thickness dependence of the 
coercive field {He oc e). As it is shown in the following 
sections, these general considerations are confirmed by 
numerical simulations. 



IV. MONTECARLO SIMULATIONS 
A. The model 

To investigate dynamic hysteresis, we perform Mon- 
tecarlo simulations based on the energy terms derived 
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above. We consider a wall of length L in a sample of 
finite dimensions. Since we are interested in macroscopic 
effects, we discretize the wall defining a minimal segment 
length a and map the zigzag wall into a particle model. 
We define the minimal elements with negative slope as 
a particle and those with positive slope as a void, as 
sketched in Fig|3| Thus we reduce the two-dimensional 
problem of a zigzag wall into a one-dimensional particle 
model evolving under the appropriate dynamic rules. 
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FIG. 3; The mapping. Two examples of zigzag configuration. 

The Montecarlo dynamics is implemented by choosing 
randomly an active pair of nearest neighbor elements, 
i.e. a particle-void or a void-particle pair, and trying to 
exchange their positions between each other. This rule 
corresponds to allow only the motion of segments with 
down-up or up-down slope, preserving the zigzag (solid- 
on-solid) structure of the wall. Once a possible displace- 
ment has been attempted, we calculate the energy dif- 
ference (see Eq. between the starting configuration 
and the new one. If AE < we accept directly the 
move, otherwise the move is accepted with probability 
P = exp[-AE/kBT]. 

In practice, the various contributions to AE are eval- 
uated as follows: 

• AE„i is obtained from Eq. ||2J). To simplify the 
expression, we can use Eij — 1/rij if the parti- 
cles are not nearest neighbor, and the whole ex- 
pression in Eq. H14|l otherwise. As it is discussed 
in the appendix, Eij deviates significantly from 
1/rij only if i and j are nearest-neighbor. We 
can then perform another simplification, by absorb- 
ing the deviation from 1/rij into the anisotropy 
(nearest-neighbor) term. In summary we can set 
Eij equal to SfioAl'^e^ /rij V«,j, and renormalize 
the anisotropy term by an appropriate magneto- 
static constant. Notice that AEm is an attractive 
long range term, tending to aggregate all the par- 
ticles. This configuration would correspond to a 
single period of a zigzag with p = L/2. 

• AEan is a nearest-neighbor repulsive term that fa- 
vors configurations where particles are followed by 
voids. This prevents the formation of a zigzag with 
a wide amplitude. In practice this term is treated 
as a positive contribution if the left (right) nearest 
neighbor of the segment i {i + 1) to be flipped has 
opposite slope with respect to i (i 4- 1), and a van- 
ishing contribution otherwise. In this way we treat 
the rotation of the magnetic moments close to the 
zigzag wall, as uniformly distributed in a "band" 
surrounding the wall^. 



• AEdis represents the contribution from structural 
disorder, whose time-independent intensity is ran- 
domly extracted with Gaussian distribution with 
zero mean for every site. 

• AEext, the interaction with the external magnetic 
field Hf^xt^ is calculated from Eg. IIUI i.e. AEext = 
—2fioHextAM , where AM is the magnetization dif- 
ference between the two configurations. 

If the move is accepted we update the configuration 
and continue the process for a time interval At. In the 
spirit of the Montecarlo method, each attempt corre- 
sponds to a time step r = 1/Nact, where Nact is the 
number of active particles (i.e. those that are followed or 
preceded by a void). After each time interval At is ex- 
pired we increase the external field by AH^xt and restart 
the updating process. This rule corresponds to a field 
rate H = AHext/ At oc w. We begin the simulation from 
the M = at H^xt = state and drive the sample to 
positive and then to negative saturation. 

B. Results 

1. Frequency dependence 
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FIG. 4: (Color online) Some T = hysteresis cycles for dif- 
ferent external field time rate. 

The first issue to be addressed in dynamic hysteresis 
is clearly the effect of the external field frequency on the 
hysteresis. In Fig. 01 we show hysteresis loops obtained at 
T — Q for various frequencies. As expected from exper- 
iments and general considerations, small (high) frequen- 
cies correspond to narrow (large) cycles. To quantify this 
observation we can focus on the coercive field behavior. 
In Fig. 13 we show the dependence of He on the field 
frequency dH/dt = H. 

At T = (for a discussion of the behavior at non van- 
ishing temperatures, see the following subsection). He 
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locity V cx dM/dt follows 

v^C{H- Hpf, 

for H slightly larger than Hp, while v vanishes for H < 
Hp. Using Eq. H13|l as a constitutive law, one can readily 
show that the dynamic coercive field scales as in Eq. H12|) 
with a = l/(/?+l). The limit a = 1 corresponds to /3 ^ 
or else to a gap in the domain wall velocity as a function 
of the external magnetic field. A very sharp dependence 
of the velocity on the field is indeed observed in our model 
(see Fig. [HJ around the coercive field, where the segments 
of the zigzag wall begin to move. This is probably a 
strong pinning effect: due to the zigzag structure of the 
wall the system is trapped by strong anisotropy barriers 
and collective effects, typically leading to a continuous 
depinning, are suppressed. 
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FIG. 5: (Color online) Coercive field vs external field fre- 
quency at various temperatures, in a linear (a) and semilog 
plot (b). Every value is mediated over 1000 realizations cor- 
responding to different disorder configurations The solid lines 
on the T — curves are a linear fit. 
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FIG. 6: (Color online) Velocity of the domain wall, normalized 
to the number of unit segments of the zigzag, as a function of 
an external constant field. Each point is an average over 2000 
different disorder configurations. 



shows an increasing linear dependence on H, of the form 
He = Hp + AH, where A is a constant. This means that 
in the adiabatic limit (low frequencies) He goes to a non- 
vanishing value Hp (as clearly shown in the semilogarith- 
mic plot of FiglSti), that we can interpret as the pinning 
dominated quasistatic component due to structural dis- 
order and anisotropy, while the linear behavior of He in 
the high frequencies regime represents the domain wall 
dominated dynamic contribution. 

This result is a particular case (a = 1) of the law 

He = Hp + AH'', (12) 

suggested on an experimental basis in Ref. (but with 
a different exponent a) an d by the theory presented in 
Ref.l20l. In the model of Ref. l20l the exponent a is related 
to the scaling exponent f3 associated with the depinning 
transition of the domain wall. In particular, it is assumed 
that under a constant applied field the domain wall ve- 



2. Temperature dependence 

Another interesting issue to analyze is the effect of 
temperature T on dynamic hysteresis. Even if thermal 
effects do not seem to be relevant for most ferromagnetic 
thin films where fc^T « AE, as it is easily checked, 
e.g., comparing ksT at room temperature, where the ex- 
periments are typically performed, with an estimate of 
AE ~ Ean which could be obtained with the parameters 
of Fe/GaAs{001) given in Sec. IIIII Thermal activation 
could play a role at very low frequencies or in ultrathin 
ferromagnetic films. In general terms, the increase of the 
temperature (at reasonably low frequency, see Fig. ^ 
acts on the hysteresis cycles shape in an similar way as 
the decrease of the frequency (see Fig. [Tj). Since a tem- 
perature increase enhances the probability for the wall to 
overcome energy barriers, at high (low) temperature hys- 
teresis loops will be large (narrow) . However at very high 
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frequencies, the system is not able to readily respond to 
the external field, and the decreasing dependence of 
on increasing temperatures is violated, as shown in fig. 
|S1 This explains the crossover between the curves with 
T = 0.4 and T = 0.8 (Fig. E)). 

In Fig. [Slwe show the dependence of He vs H at vari- 
ous non-zero temperatures. It is interesting to note that 
our simulations predict that the high T behavior of He 
vs H deviates from the linear behavior established at 
T = 0. This can be understood from general consider- 
ations, since when T > thermal activation will lead 
creep domain wall motion even for H < Hp. Hence, ac- 
cording to the theoretical analysis presented in Ref. |^ 
at low frequencies the dynamic coercive field will deviate 
from Eq. ((T^ and decay as He ~ l/[log(ij)]^/'^ where fi 
is a creep exponent. This result is consistent with Fig. 
|5Jd, although the limited scaling range does not allow for 
a quantitative confirmation. 
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FIG. 8: (Color online) Coercive field vs temperature for var- 
ious external field time rates. Every value is mediated over 
1000 realizations corresponding to different disorder configu- 
rations. 
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FIG. 7: (Color online) Some liysteresis cycles for various tem- 
perature. Tfie external fixed field time rate is At = 100. 



3. Thickness dependence 

Finally, we address to the film thickness role in dy- 
namic hysteresis. As it can be seen in Eq. ^2 for suf- 
ficiently large e above the purely disorder dominated 
regime discussed in section IIII Bl the energy barrier in- 
creases linearly with e at T = 0, and the coercive field 
does so as well. 

Turning to a more quantitatively discussion, at T = 
the results of Eg. lTTI are confirmed by the simulations for 
various e summarized in Fig. ^ In the two upper curves, 
where the disorder energy term is not negligible, we no- 
tice three regimes: a divergence proportional to at 
very low e, due to the disorder term and linear regime 
for high e values, due to the dipolar term (the regime 
independent from e that would be due to the anisotropy 
energy and the external field, cannot be seen clearly in 




FIG. 9: (Color online) Coercive field vs film thickness for 
various values of the a of the disorder Gaussian distribution. 
Every value is mediated over 200 realizations corresponding 
to different disorder configurations. 

that figure). Moreover, from the lower curve of Fig. 
(the one without disorder) we confirm that the low-e di- 
vergence is due to the disorder term. Otherwise, the 
thickness dependence of the coercive field does not affect 
the linearity of the frequency dependence of He, which 
remains valid for every e value, as is seen in Fig. 1101 (at 
T = 0). 



V. COMPARISON WITH EXPERIMENTS 

The linear dependence of He vs H is consistent with 
the experimentsSc, at least in the high frequency regime, 
as can be seen from Fig. even though there the au- 
thors try to identify two dynamical regimes in a semi- 
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FIG. 10: (Color online) Coercive field vs frequency H a.t T = 
for two different thicknesses. The linear behavior is not 
affected by the thickness. Compare the present figure with 
the experimental results reported in Fig. 1111 



show a deviation from the linear dependence. Even if this 
deviation is quite small and could be due simply to the 
experimental error, our results suggest that it could be 
interpreted by means of thermal activated motion. This 
interpretation seems to be supported by the fact that the 
deviation is stronger for the thinner sample, as it would 
be expected if thermal effects where the cause. We notice 
as well that the data reported in Ref.^.3. for even thinner 
films of Fe/GaAs show a logarithmic decrease of He at 
low frequency, not observed for thicker samples where 
He tends instead to a constant as in Fig. Finally 
our model accounts as well for the fact that the value 
of He for the thinner sample (in our notation, the one 
with e = 55A) is smaller than the value measured for the 
thicker one (e = 250A) at the same field frequency H. 
As a word of caution it is important to remark here that 
the only sample for which a zigzag domain structure is 
clearly reported is the one with 250A thickness. 



VI. CONCLUSIONS 



logarithmic plot, attributed to two different dynamic pro- 
cesses, namely domain nucleation and domain wall prop- 
agation. This is probably an artifact of the logarithmic 
scale. Our conclusion is rather different: as it is quite 
clear from figure lTTl the experimental results-^- on the fre- 
quency dependence of the coercive field could be properly 
explained by a model based on a single propagating do- 
main wall, which is the only dynamic process that we 
have considered, without invoking domain nucleation. 




□ E=5.5nm - 

ftoi 0.1 1 10 100 1000 

H (kOe/s) 

FIG. 11: (Color online) Coercive field as a function of the 
external field frequency from experiments on Fe/GaAs thin 
films for two different thicknesses. The data is fitted properly 
by linear functions (solid lines), but when plotted in semi- 
logarithmic scale they may suggest a crossover at ~ 20. 
In the inset: a linear plot of the frequencies up to lOOfeOe/s. 
The deviation from linear behavior is quite small and could 
be do to thermal activation. Data are obtained from Ref. ^ 

In the low frequency regime, the experimental data 



The dynamics of two dimensional ferromagnetic sys- 
tems is still under debate, both on the theoretical and 
the experimental side. A crucial issue is the description 
of the dynamic hysteresis, that is related to power losses 
and thus plays an important role in several technological 
applications. Here, we have analyzed ferromagnetic thin 
films with zigzag domain walls, arising when the mag- 
netization vectors in two nearest-neighbor domains meet 
head-on at the wall that separates them. To investigate 
dynamic hysteresis, we have studied the motion of zigzag 
domain walls by developing a simple discrete model based 
on the interplay between dipolar and anisotropy energy 
contributions, in presence of structural disorder. Un- 
der some simple approximations one can estimate some 
experimentally relevant quantities, such as the typical 
zigzag half-period and the coercive field, which turn out 
to be in quantitative agreement with experimental obser- 
vations. 

Although quite simplified, our model allows to recover 
the behavior of coercive field He in dynamic hysteresis. 
We have studied the dependence of He on the applied 
magnetic field frequency H at T = Q and found that the 
coercive field scales as He = Hp + AH. This linear be- 
havior is in good agreement with experiments^, which we 
can thus explain by means of pure domain wall propa- 
gation model, without the need to invoke other dynamic 
processes as domain nucleation. We have also simulated 
hysteresis at T > 0, even if this case is probably not rel- 
evant for most ferromagnetic thin films, where thermal 
effects are negligible. We show that high temperature at 
low frequency induces narrow loops and the coercive field 
decreases with respect to the T — case, while at high 
frequency the situation is less intuitive due to the delay 
between the system response and the external driving 
field. We have also studied the dependence of the coerci- 
tivity from the film thickness e. The behavior indicated 
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by the simulations is explained by simple analytical con- 
siderations. For small disorder, we find that the thickness 
does not affect the rate dependence of He &t T = 0. It 
is interesting to remark that our model could be appli- 
cable to ferroelectric materials which are known to show 
as well zigzag domain walls. It would be very interesting 
to compare our results with experiments in this case as 
well. 
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VII. APPENDIX 



and 



Analitical calculation of magnetostatic energy: the re- 
sult of equation |2Jl is given by 



8M^e^ X 

{6{mi, rrij ) [gi {rm , - qj , nij , jp, ip) 
-gi{mi,qi - qj,mj,jp, {i + l)p) 
-gi{m,,q, - qj,mj, {j + l)p, ip) 
+gi{mi,q, - qj,mj, (j + l)p, (i + l)p)] 
+6{mi,-mj) [g2{nii,qi - qj,mj,jp,ip) 
-92{mi,qi - qj,mj,jp, {i + l)p) 
-92{mi, q^ - qj,mj, {j + l)p, ip) 
+92{mi, - qj,mj, (j + l)p, {i + l)p)]} , 



(14) 



where 



gi{m,q~ q',m,x,x') 



1 



n a(r-r') 
r - r In 



1 + rn^ 

r — r'l — a(r — r') 



gi{m,q- q',m,x,x') 



r — r'l + a(r — r') 

V q^q' 
{—X + [x — x') \n{x — x')) 



if q = q' 



52 ('Ti, q — q' , m, x, x') — b(r — r') In 



r — r + a r — r 



+6(r' - r) In 



r — r'l — a(r — r') 

|r - r'l + c(r - r') 
|r — r'l — c(r — r') 



a(r — r') 
&(r-r') : 
c(r — r') : 



x' — X + 'm?[x' — —x) — m'{q — q') 



Vl + m2 



q — q' + 2mx 



4to\/1 + rn^ 



X — x' + m?'{x + x') + m{q — q') 



1 -f m? 



and d{m,m') is Kronecker's delta. 

In Fig. Elwe plot the function Eij (for unitary 8M^e^) 
as a function of the distance between the centers of mass 
of the segments i and j. As could be seen, the value 
of magnetostatic interaction energy is mainly propor- 
tional to l/|r — r'l = 1/rij for each pair of i and j but 
the nearest-neighbour, where magnetic charges are suffi- 
ciently close to each other to experience the very shape 
of their spatial distribution. 
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FIG. 12: (Color online) Interaction magnetostatic energy be- 
tween two generic segments with parallel or opposite slopes 
as a function of the distance between their centers of mass. 



